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STABILITY OF A REGULARIZED NEWTON METHOD WITH TWO POTENTIALS 
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Abstract. In a Hilbert space setting, we study the stability properties of the regularized continuous Newton 
method with two potentials, which aims at solving inclusions governed by structured monotone operators. 
The Levenberg-Marquardt regularization term acts in an open loop way. As a byproduct of our study, we 
can take the regularization coefficient of bounded variation. These stability results are directly related to 
the study of numerical algorithms that combine forward-backward and Newton’s methods. 


1. Introduction 

Throughout this paper iJ is a real Hilbert space with scalar product (•,•) and norm ||a;|| = yj (a;, x) for 
X £ H. We make the following standing assumptions: ip and ip are two functions which act on H and satisfy 

• ip : H —> R U {+oo} is convex lower semicontinuous, and proper; 

• Ip : H —R. is convex differentiable, and Vip is Lipschitz continuous on the bounded subsets of H. 

• p + Ip is bounded from below on H. 

We are concerned with the study of Newton-like continuous and discrete dynamics attached to solving the 
the structured minimization problem 

(T) min{(^(a;)-I-'0(a:) : x G H} . 

Note the asymmetry between ip, which may be nonsmooth, with extended real values, and ip which is 
continuously differentiable, whence the structured property of the above problem. Indeed, we wish to design 
continuous and discrete dynamics which exploit this particular structure and involve p via implicit operations 
(like resolvent or proximal operators) and ip via explicit operations (typically gradient-like methods). So 
doing we expect obtening new forward-backward splitting methods involving Riemannian metric aspects, 
and which are close to the Newton method. This approach has been delineated in a series of recent papers, 
m, i, i, In this paper we are concerned with the stability properties with respect to the data (A, xq, 
Vo, ...) of the strong solutions of the differential inclusion 

v(t) G dp{x(t)) (la) 

X{t)x{t) -I- v{t) -I- v{t) + 'Vip{x{t)) = 0 (lb) 

x(0) = a;o,u(0) = uo- (Ic) 

Let us now make our standing assumption on function A(-): 

(2) A : [0, - 1-00 [—]0, - 1 - 00 [ is continuous, 

(3) and absolutely continuous on each interval [0, 6], 0 < b < -|-oo. 

Hence X{t) exists for almost every t > 0, and A(-) is Lebesgue integrable on each bounded interval [0,6]. We 
stress the fact that we assume X{t) > 0, for any t > 0. By continuity of A(-), this implies that, for any 6 > 0, 
there exists some positive finite lower and upper bounds for A(-) on [0,6], i.e., for any t G [0,6] 

(4) 0 < Xly^rain ^ A(t) ^ ^ -t-OO. 

Our main interest is to allow X{t) to go to zero as 6 —> -l-oo. This makes the corresponding Levenberg- 
Marquardt regularization method asymptotically close to the Newton’s method. 

Let us summarize the results obtained in [5]. Under the above assumptions, for any Cauchy data xq G domclt/J 
and Vo G dp{xo), there exists a unique strong global solution {x (•), v (•)) : [0, -|-oo [—H x H oi the Cauchy 
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problem (ITal) -(fTH. Assuming that the solution set is nonempty, if X{t) tends to zero not too fast, as 
t —+ 00 , then v (t) —> 0 strongly, and x (t) converges weakly to some equilibrium which is a solution of 
the minimization problem (CP). By Minty representation of dtp, the solution pair {x (•), v (•)) of (ITall - dT^ can 
be represented as follows: set fi{t) = then for any t G [0, +oo), 

(5) a;(t) = prox^(t)^ (z(t)); 

(6) r;(<) = (2:(t)) , 

where z (■) : [0, +oo[— H is the unique strong global solution of the Cauchy problem 

(7) z{t) + {p{t) - /i(t))(z(t)) + {z (t))) = 0 

(8) 2;(0) = a:o+ At(0)uo- 

Let us recall that prox^,^ is the proximal mapping associated to fip. Equivalently, prox^^ = (/ + 
is the resolvent of index /r > 0 of the maximal monotone operator dp, and is its Yosida approximation 
of index p > 0. 

Let us stress the fact that, for each t > 0, the operators prox^^^j,^ : H —> H, S/p^f^t) ■ H — H are 
everywhere defined and Lipschitz continuous, which makes this system relevant to the Cauchy-Lipschitz 
theorem in the nonautonomous case, which naturally suggests good stability results of the solution of (fTal) - 
© with respect to the data. 

This paper is organized as follows: We first establish a priori energy estimates on the trajectories. Then 
we consider the case where A is locally absolutely continuous. Note that it is important, for numerical 
reasons, to study the stability of the solution with respect to perturbations of the data, and in particular of 
A which plays a crucial role in the regularization process. In Theorem 13.II we prove the Lipschitz continuous 
dependence of the solution with respect to A. Moreover, the Lipschitz constant only depends on the norm 
of the time derivative of A. Finally, we extend our analysis to the case where A is a function with bounded 
variation (possibly involving jumps). We use a regularization by convolution method in order to reduce to 
the smooth case, and then pass to the limit in the equations. So doing, in Theorem 14.31 and Corollary 14.41 
we prove the existence and uniqueness of a strong solution for (la)-(lc), in the case where A is a function 
with bounded variation. 


2. A PRIORI ESTIMATES 

The linear space H x H is equipped with its usual Hilbertian norm ||(^, C)|| = '\f ll^ll^ + IICII^ • In this 
section we work on a fixed bounded interval [0,r], and following assumption Q, we suppose that there 
exists some positive constant cq such that 

(9) 0 < Co < X{t) for all t G [0,T]. 

We will also assume that V?/; is L^-Lipschitz continuous. Indeed, this is not a restrictive assumption since 
one can reduce the study to trajectories belonging to a fixed ball in H. We will often omit the time variable 
t and write x,v.... for x (t), v (t).... when no ambiguity arises. 

In the following two Propositions we denote by {x (•), v (•)) : [0, +oo) H x H the unique strong global 
solution of the Cauchy problem (TTall - dT^ . 

Proposition 2.1. Let {x,v) be the strong solution of system ila\) - ilc\) on [0,r], T > 0. Then 

(10) ^ ||i;(t)f < ^ ((:^ + ^/>)(xo) - inf((^ + '0)) . 

(11) ||a:||L~(o.T;R) < Ikoll + ((<P + if){xo) - inf(‘/? + V")) " ■ 

Proof, a) For almost every t > 0, x{t) and v{t) are well defined, thus 

{x{t),v{t)) = lim -^{x{t + h) - x{t),v{t + h) - v{t)). 
h^o n‘‘ 

By equation (fTal) . we have v{t) G dp{x{t)). Since dp : H ^ H is monotone 

{x{t + h) — x{t),v{t + h) — v{t)) > 0. 
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Dividing by and passing to the limit preserves the inequality, which yields 

(12) (x{t),v{t)) >0. 

By taking the inner product of both sides of m by x{t) we obtain 

A(t)||i;(t)f + {x{t),v{t)) + {v{t),x{t)) + {V'tlj{x{t)),x{t)) = 0. 

By (TT^ we infer 

(13) A(t)||i:(t)f + {v{t),x{t)) + {V'tp{x{t)),x{t)) < 0. 


Noticing that x and v are continuous on [0,T], hence bounded, and that A is bounded from below on [0,T] 
by a positive number, one easily gets from (1131) that 

(14) x(^L^(0,T]H). 

For our stability analysis, we now establish a precise estimate of the norm of x. By the classical derivation 
chain rule 

(15) = {V-il){x{t)),x{t)). 

We appeal to a similar formula which is still valid for a convex lower semicontinuous proper function ip : 
—?► RU {+oo}. Notice that 

i) v(t) € dp{x(t)) for almost every t G [0,T]; 

ii) V is continuous on [0,T], and hence belongs to L'^{0,T\ H); 

in) X G by IT^ . 

By i), ii), Hi), conditions of El Lemma 3.3] are satished, which allows to deduce that t p{x (t)) is 
absolutely continuous on [0,T], and for almost t G [0,r], 

(16) ^p{x{t)) = {v{t),x{t)). 

Combining (IT^ with (ITSl) and (fTO we obtain 


(17) A(t)||a;(t)f + ^ {p{x{t)) + V'(a;(t))) < 0. 

By integrating the above inequality from 0 to T we obtain 

(18) [ X{t)\\x{t)\\'^dt+{p + tp){x{T)) < {p + tjj){x{0)). 

Jo 

Since (p is bounded from below on and A is minorized by the positive constant cq on [0, T] (see ((H])), 
we infer 

pT ^ 

(19) ||i;(t)fdt < — (((/? +V')(a:o)-inf(<p+ V')) • 
b) Since a;(-) is absolutely continuous on bounded sets, for any t G [0,r] 


( 20 ) 


?(t) =Xo+ x{T)dT. 

Jo 


Passing to the norm, and using Cauchy-Schwarz inequality yields 


( 21 ) 


( 22 ) 


lk(t)|| < lixoll + / ||i:(T)||dr 

Jo 


< \\xo\\ + Vt ||a;(r)fdT 


Combining the above inequality with (11911 gives 


( 23 ) 


lk(0ll < Ikoll + \ — ((7> + i’)ixo) - inf(i^ + 

V Cn V H 


This being true for any t G [0, T] 


(24) 


L-(o,T;H) < Ikoll + \/ — {(<P + '0)(a:o) - inf(‘/? + 


Let us now exploit another a priori energy estimate. 

Proposition 2.2. Let {x,v) be the strong solution of system / ftali -i lTel] on [0,T], T > 0. Then 


(25) 


2T^L'^ 

\\v{t)\\^dt < ||uof + 2T||VV'(xo)f + [{t + ^){xo) - inf(</5 + ^/>)J 


Co 

(26) ||c||L“(o,T;ff) < Ikoll + V^\\yif{xo)\\ + J—TL^ + iP){xq) - inf((^ + 

V Co V n 

Proof. By taking the scalar product of (fTb| by v{t) we obtain 

X{t){x{t),v{t)) + ||u(t)f + {v{t),v{t)) + {yif{x{t)),v{t)) = 0. 

By uni) we infer 

(27) 

Hence 

(28) 


||u(t)f + {v{t),v{t)) + {V'ip{x{t)),v{t)) < 0. 


^ l|V^(a;(i))||||c(0ll 


2 dt 


(29) 

which implies 

(30) 


+ l\mr 


^IKi)f <||VV-(a;(t))f. 


By integrating the above inequality we deduce that, for any t G [0,T] 


(31) 


f ||u(r)f dr + ||u(t)f < lluof + T || VV^(x(r))f dr. 
Jo Jo 


Since is L^-Lipschitz continuous 

(32) ||VV’(x(r))|| < ||VV'(a:o)|| + Li^W^ir) - xo||. 

A careful look at the proof of (ETl) gives the more precise estimate 

Ft 


(33) Ik — a;o||L°o(o,T;i/) < 
Hence for any t G [0, T] 

(34) ||VV-(x(r))|| < ||V^(a:o)||+L^ 
Combining (I5T]) with (IMl) gives 


ip + 'ip){xo) - mi{(p + 

iri 


(/? + V’)(a:o) - inkT’■ 

H 


(35) ||u(r)f dr + ||u(t)f < ||uof + 2T||VV’(a:o)f + {(P + ^)ixo) - + i/')) . 


As a consequence 

(36) f 

Jo 

and 


\\v{t)fdt < ||uok + 2T||VV'(a:o)k 4-((7> + '<P){xo) - inf(v3 + '0)) 


(37) lkllL“(0,T;i/) < Ikoll + v^||V0(xo)|| + W— ({(fi + tj;){xo) -inf((/j + 

V Co V u 

which ends the proof. 


□ 


□ 


We can now deduce from the two preceding propositions an a priori bound on the L°° norm of x and v. 
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Proposition 2.3. 


(38) l|i||L“(0,T;ff) <-1- 


CO 




Co 


ip + tp){xo) - inf(iy9 + 


H 


(39) 


\i\\L°°{0,T-,H) < 


(1 + v^)||VV^(a;o)|| + (^'^ +^)L^ (^{ip + 'ijj){xo) - ini{p + tp)^ " 


Proof, a) Let us return to the equation obtained by taking the inner product of both sides of (HB by i{t) 
A(<)||a;(t)f + {x{t),v{t)) + {v{t),x{t)) + {\7'ilj{x{t)),x{t)) = 0. 

By (1121) . and A is minorized by the positive constant cq on [0,T],we infer 


co||a:( 


{v{t) + V'ip{x{t)),x{t)) < 0. 


Hence 


||a;(t)|| < —||u(t) + VV'(a;(t))|| 

Co 

< —(ll'c(<)ll + II V?/'(a:o)|| + L^\\x{t) - a:o||). 
Co 


l^c||L~(0,r;ff) < 


By (1551) and (1571) we deduce that 
1 

Co 
1 

Co. Co V Cq 

lluoll 1 + v^ 


||uo|| + v^||V^(a:o)|| + \ —TL^p ({p + i’){xo) - inf((^ + 


Co 


H 


< 


— ||VV'(a;o)|| + — ((V5 + '0)(a;o) - inf((^ + 

Cn Cn V Cn V H 

{V2T + VT)L^ / 


l|VV'(xo)|| + 


1 ((/5 +'0)(a;o) - inf((/3 


Co Co 

b) Let us return to the equation obtained by taking the inner product of both sides of ([ib| by v{t) 
X{t){x{t),v{t)) + ||u(t)f + {v{t),v{t)) + {Vilj{x{t)),v{t)) = 0. 

A similar argument as above yields 

||f'(i)|| < ||u(t) + VV'(a;(t))|| 

from which we deduce the result. 


□ 


Let us enunciate some straight consequences of Proposition [23] • 

Corollary 2.4. The following properties hold: for any 0 < T < +oo 
(1) t I—> x{t) is Lipschitz continuous on [0,r] with constant 


Co 




Co 


p + if){xo) - ini{p ■ 
H 


(2) 1 1 — v{t) is Lipschitz continuous on [0,r], with constant 

{V2T + VT)L^ / 


+ (l + v^)||VV'(xo)|| + 


1 ((/? +'0)(a;o) - inf(v3 


3. Stability results 

In the next theorem we analyze the Lipschitz continuous dependence of the solution (x, v) of the Cauchy 
problem (llal) - (|lcl) with respect to the function A and the initial point (xo,vo). Our demonstration is based 
on that followed in 0 Theorem 3.1]. 
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Theorem 3.1. Suppose that X,r] : [0,T] —> [co,+c»[ are absolutely continuous functions, with T > 0 and 
Co > 0. Let (x,v), (y,w) : [0,r] — > H x H be the respective strong solutions of the inclusions 

(40) Xi + i) + v + Vtp (x) = 0, V € dp (x), x (0) = xg, v (0) = vq, 

(41) r]y + w + w + X/'ip{y) = 0, w € dp{y), y (0) = yo, w (0) = wq. 

Define 9 : [0, T] —> K, for each t G [0, T\, by 

e{t) = sj cl ||a;(t) - y{t)\f + ||w(<) - w{t)f. 

Then 

(42) 

ll^'l|L“([o.r]) < 

with 


A( 0 ) + 77(0) 11 11 11 .. ^ w \ II 

-2- - 11^0 - 2 /o|| + ||uo -u'oll + ^l|A- ? 7 lUi([o.T]) 


X exp 




2cn 


+ T (1 + ^) 
Co , 


^ ^ IKII + lkoll ^ + IIVV^(yo)||) + 


Co Co 

{V2T + VT)L^ 


p + 'if){xQ) - ini{p ■ 

M 


+ ({T + i’)iyo) - inf(((? + 


In particular, if xg = yg, vg = wg, then 


(43) 

with 


C 

||6'||l~([0,t]) < -^l|A-?7 ||li([0.t]) X exp 


11 ^ + ^IIli([o,t]) 


2 Co 


+ T (1 + ^^) 
Co , 




Cg Cg 

Proof. To simplify the exposition, define 7 : [0,T] 


p + ii){xg) - inf((p ■ 

H 


X + T] 


7 = 

Using the assumptions X,r/> cg and the monotonicity of dp we conclude that for any t G [0,r] 

Cg < 1 , {x - y,v - w) > 0 . 

Therefore 

(44) 9 <\J-i'^\\x- y\f + ||u - w||^ < II 7 {x - y) + {v - u;)|| . 

Let us show that 9 satisfies a Gronwall inequality. Let us start from 

d 


dt 


[7 {x - y) + {v - w)] = j {x - y) + j {x - y) + {v - w). 


In view of (|40)) and dUD, 

{v — w) = —Xx + rjy — (v — w) — (X/ip (x) — V'0 (y)). 

Combining the two above relations gives 

{x - y) + {v - w)] = j {x - y) + ^ {x - y) - Xx + py - {v - w) - {Vfj (x) - Vip (y)), 

= iix-y) + (x + y)- {v-w)- {ytp (x) - Vf) (y)). 

Since 7 , x, y, w, v are absolutely continuous, the function ^ {x — y) + {v — w) is also absolutely continuous. 
As a consequence, integration of the above inequality on [0, s], for s G [0,T], yields 

[7 {x-y) + {v- w)] (s) -b{x-y) + {v- w)] ( 0 ) = 

iv->^){i + y)dt+ f bl{x-y)-{v-w))dt- [ {X/f; (x) - Vip (y)) dt. 

^ Jo Jo Jo 
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Passing to the norm 

II [7 (x - y) + (z; - w)] (s) II < 7 (0) ||a;o - ?/o II + Iko - Wo II 

(45) f \^-v\\\i + y\\dt+ [ \\'^{x-y)-{v-w)\\dt+ [ WWip (x) - Wip {y) \\dt. 

^ Jo Jo Jo 

By Lipschitz continuity property of V^, and definition of 9^ we have 


(46) 

On the other hand 


||Vz/> (x(t)) - Vz/> (y(t)) II < L^\\x{t) - y{t)\\ 

< —e(t). 

Co 


Il7(a; - y) - (^^-w)|| = |7l^ Ik - y|| +|k-w|| +2{co{x-y), — {w-v) 

Co 

l7l^ 


< I 7 I lk-y|| +||z;-w|| +Co|k-y|| + — ||z; - zi;|| 


(47) 


iT-r 


+ 1 r 


Combining (HSl) . (1451) . (ITTl) and 9 < ||7 [x — y) + {v — w)||, we obtain 


6 »(s) < 7 (0) Iko - yoll + Iko - wo|| + ^ / |A - yl ||x + y|| dt+ f 9{t) W ^ + 1 + 

2 7o Jo V ^0 Co 


< 7 (0) Iko - yoll + Iko - woll + 2 lk “ ylki([o,T])lk + y||L“([o,T]) 

Applying Gronwall’s inequality yields 


9{t) 


I 7 I 

—+ 1 + 

Co Co 


dt 

-7/y> 


dt. 


9{s) < 


7 (0) Iko - yoll + Iko - Woll + -||A - y||Li([ 0 .T])lk + ylk“([ 0 .r]) 


X exp 


hl + i + A 


Co 


Co 


dt. 


Combining this estimation with the bound on the L°° norm of x and y (see Proposition (12.31) ') gives 


lkl|L<-([0,T]) < 


c, 


7 (0) Iko - yoll + Iko - Woll + 17IIA - y||Li([o,T]) 


X exp 


hl + i + ta 


Co 


Co 


dt 


with 


^ ^ Ikoll + llwoll ^ ^ IIvV.(yo)||)+ 

Co Co 

{V2T + VT)L^ 


H 


if + 'tp){xo) - inUif + + (<^ + '^){yo) - inf((^ + 


H 


Equivalently 
lkl|L‘-([0,T]) < 


A( 0 ) + 77(0) II II II II ^11 \ II 
- ^- Iko - yoll + Iko - woll + —||A - y|ki([o,T]) 


xexp 


l|A + ylkqio.T]) 


2 co 


+ T (1 + ^) 

Co 


□ 


Remark 3.2. In the cas tp = 0 we recover exactly the same stability results as [i Theorem 3.1]. 

Remark 3.3. It is worth noticing that, besides the Lipschitz continuous dependence with respect to A of 
the solution (x, v) of (fTal)-(fTB. Theorem (ED also provides its continuous dependence with respect to the 
initial data (xo,uo)- More precisely, if (x„,z;„) (resp. {x,v) ) is the solution of (fTa1) - (IT^ corresponding to 
the Cauchy data (xon, x’on)(i'esp. {xo,vo)), as a direct consequence of (H^ we obtain for all T > 0 

{dip 9 (xon,von) — (xo,z;o)) (xn,Vn) — t (x,v) Uniformly on [0,T]. 

Note also that, by contrast with semigroup generated by a subdifferential convex there is no regular¬ 
izing effect on the initial condition: there is no way to define a solution of (fTa 1 ) - dT^ for xq G domcit/?\dom 9(/3 
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since in that case, for any approximation sequence xon G doni9(/3 and von € d(p{xon), one has lim||r;o„|l = 

n 

+ 00 , which would imply the blow-up of the sequence {xn,Vn) as n —> -foo, on any finite interval. 

Remark 3.4. Another approach is to study the equivalent problem ©-([S]), based on the known stability 
results for the Cauchy-Lipschitz problem. Although conceptually simple, this approach seems more technical. 

4 . Bounded Variation Regularization Coefficient A(-) 

Let us suppose that A(-) : [0,r] —]0, oo[ is of bounded variation on [0,T], where T > 0. That is 
TV (A, [0,T]) < - 1 - 00 , where TV (A, [0,T]) is the total variation of A on [0,T] : 

p 

TV (A, [0,T]) = sup^ |A(Ti) - A(Ti_i)|, 

i=l 

the supremum being taken over all p G N and all strictly increasing sequences tq < ti < • • • < Tp of points 
of [0,T]. Function A may involve jumps. We also suppose that A is bounded away from 0: 

inf A([0,T]) > 0. 

The following lemmas which are proved in [3 Lemma 3.1, Lemma 3.2] gather some classical facts concerning 
the approximation of functions of bounded variation by smooth functions together with some technical results 
useful for sequel. 

Lemma 4.1. Let A : [0,T] —> ]0, oo[ be of bounded variation on [0,T]. Then there exists a sequenee 
(An)„gi^, with Xn G C°° ([0,T]) for each n G N, such that 

i) inf A([0,r]) < A„ (t) < supA([0,r]), Vt G [0,T], Vn G N. 

In particular, A„ > 0 */ A > 0; 

ii) Xn —> X in (0, T) for any 1 < p < oo; 

lii) rV(A„, [0,T]) = fj' |A„ (t) Idt < TV(X, [0,T]). 


Lemma 4.2. Let Zn,z G C ([0,T] ,H) be such that Zn —> z uniformly and {zn)n is L-Lipschitz continuous 
for some positive constant L independent o/n G N. Let Xn —>■ A in (0,T). Then XnZn converges weakly 
to Xz in L^ (0, T; H). 


We can now state the main result of this section. 


Theorem 4.3. Let X : [0, T] — > ]0, oo[ be of bounded variation on [0, T], and suppose cq = inf A ([0, T]) > 0. 
Let xq G domdip and vq G d(p(xo), vq ^ 0. Then there is existence and uniqueness of a strong solution 
(x,v) : [0,T] — > H X H of the Cauchy problem 

(48a) V (t) G dp {x (t)), 0<t<T 

(48b) X{t) X {t) + V {t) + V {t)+ V'ijj {x (t)) = 0, a.e. 0<t<T, 

(48c) x{0) = xo, v{0)=vo. 


Proof Existence: According to Lemma im there exists a sequence of functions (A„)„gj^ in (7°° ([0,T]), 
with A„(t) > Co, which converges to A in L^ (0,T) for p > 1, and satisfies the condition 

rT 


(49) 


A„ (t) dt < TV (A). 


For each n G N there exists a unique {x„,Vn) solution of the differential inclusion 


(50a) Vn (t) e dp {xn{t)), 0<t< T, 

(50b) Xnit)Xnit) +Vn{t) +Vn{t) +V->p {Xnit)) = 0, O.C. 0 < t < T, 

(50c) XniO)=Xo, VniO)=Vo- 


We will show that (a;„,u„) converges uniformly to a solution of (l48aM8cp . 





Consider, in C ([0, T], iJ x H), the norm 


I = max 
“ ie[0,T] 


\Jcl \\z{t)f + ||w(t)||' 


It is equivalent to the sup norm in C {[0,T] ,H x H). Using Theorem 13.11 (l43l) and (l4^ we deduce the 
existence of a constant C (which is independant of n) such that for any n, m we have 


i,u„) - {Xm,Vm)\\c„ < C exp 




2 Co 


+ T(1 + I ||A„ - 


< Cexp 


^+ ^(1 + 11'^" ~ '^"i|Ili([q_t]) ■ 


Since the sequence (A„),^gj^ converges to A in (0,r), we deduce that is a Cauchy sequence with 

respect to the sup norm. Therefore, {xn, Vn) converges uniformly to some continuous (x, v) : [0, T] — HxH. 
Since is continuous we also obtain 


V^(x„) —>■ V^(x) uniformly on [0,T]. 

Moreover, Lemmaensures that A„i;„ and converge weakly to Xx and i) in (0,T; H). Hence letting 
n —>• +0O in (|50bl) gives 

Xx + V + v + yip (x) = 0, 

that’s (I48bl) . Finally, since the graph of d(p is closed we obtain, for all t G [0,T], v (t) G d<p{x{t)) that’s 
(I48al) . Hence {x,v) is a solution of (l48al) - (l48cD . 

Uniqueness: We adapt the proof of [5], using differential and integral calculus for BV functions which 
involves differential measures. 

Define A“ : [0,T] i —> [co,+oo[ by 


A- (0) = A(0), 

0 < t < T : X~ (t) = lim X{t — e). 

e>0,e^0 

Let (x,v ), {y,w) : [0,T] ^ H x H he two strong solutions of (I48al) - (l48bl) - (l48c|) . Explicitly 

Xx + v -\-v + yip {x) = 0 a.e.; v (t) G dp {x (t)) Vt; x (0) = xq, v (0) = vq, 
Xy + w + w + yip {y) = 0 a.e.; w (t) G dp {y {t)) Vt; y (0) = Xq, w (0) = vo 


Since A = A a.e., we also have 

X~x + ii + V + yip (x) = 0 and X~y + it + w + yip (y) = 0 a.e. 


and consequently 


(51) X~ {x — y) + {v — it) + V — w + yip (x) — yip (y) = 0 a.e. 

In terms of differential measures on [0,T] we have [221 Proposition 11.1] 

(52) d [A“ {x — y) + {v — ic)] = X~d {x — y) + {x — y) dX~ -\- d{v — w) . 

Integrating the left hand term on [0, s[ and taking the initial condition into account, we obtain for s G [0, T\ 
[m Corollary 8.2] 


f d[A {x-y) + {v-w)]=X (s)(x(s) 
J[oA 


y{s)) + (u(s) - w (s)). 


(53) 





















Now integrating the right hand term of (15^ on [0,s) and taking (ISTl) into account, we get 
/ [A“(i (x — y) + (x — y) dX~ + d{v — rc)] 

d[0.s[ 

= / [X~d {x — y) + d {v — w)] + / {x — y)dX~ 


/ [a {x — y) + i) — w] dt + (x — y) dX 

^\0,s[ '^[0,s[ 

/ [v — w + X/tjj {x) — Vil: (y)] dt+ {x — y) dX~. 

Mo.sl J\0,s\ 


(54) 

From (I52I53I54I) we deduce 

A“ (s) (x (s) — y (s)) + {v (s) — w (s)) = / {x — y) dX~ — [u — w + V'0 (x) — (y)] dt. 

^fO,s[ ‘^fO,s[ 


Whence 
(55) 


A (s)(x(s)-?/(s)) + (u(s)-u;(s))|| < / ||x- 2 /|||dA \+ ( \\v-w\\dt+[ \\Vijj {x) - X/{y)\\ dt 

J[0,s[ '^[0,s[ J[0,s[ 

1/2 

6 (s) = ^Cq ||x (s) — y (s)|/ + Hu (s) — w (s)|/^ . The same reasoning as in Theorem 13.II yields 


Define 

0is) < IIA" (s) (x(s) - 2/(s)) + (u(s) - ■u;(s))|| , Vs e [0,r]. 

Besides we also have cq ||x (s) — y (s)|| < 9 (s), ||u (s) — w (s)|| < 0 (s), and we have that X/ijj is the gradient 
of a convex, continuously differentiable function ip : H —> M, and by Lipschitz continuity property of V'i/;, 
and definition of 0, we have 

Wx/ip (x(<)) - X/ip {y(t)) II < L^||x(t) - y{t)\\ 


(56) 

Hence, with (l55l) 

(57) 


< ^9{t). 

Co 


0 (s) < — f 9\dX 1+ f 9dt-\ - [ 9dt = f 

Co A[o.s[ JlO,s[ Co J[o,s[ 0[c 

where dy, denotes the nonnegative measure MA“| + ^1 + dt. 


9dy, 


If 0 ^ 0 on [0, T] , define to = inf {t G [0, T] ,9 (t) > 0}. Note to < T and 0 (to) = 0, since 0 is continuous. 
With (1571) we then have 

(58) 0 (s) < [ 9dy, to < s < T. 

•']to,s[ 

In view of dy = 0 and of the right continuity at to of t —>■ dy, [22l Proposition 9.1 ] there exists 

some ti G ]to,T] such that dy < 1/2. Let M be an upper bound of 0 on [0,ti]; from (l58)) we deduce, 

for s € ]to,ti] 

f r M 

9 (s) < M dy<M dy < —. 

A]to,s[ AtoXil ^ 

Hence M/2 is also an upper bound of 0 on [0, ti], which necessarily entails M = 0 and 0 = 0 on [0, ti]. But 
this is contradiction with the definition of to- Hence 0 = 0 and (x,u) = {y,w) on [0,T]. □ 

Theorem 2.2 has a natural global version formulation; 

Corollary 4.4. Suppose that X : [0, oo[ —> ]0,oo[ is of bounded variation on [0,r] and inf A([0,r]) > 0 for 
any T < oo. Let vo G dip{xo) and vo ^ 0. Then there is existence and uniqueness of a strong solution 
(x,u) : [0,oo[ —> H X H of the Cauchy problem 

(59) Ax + u + u + X7ip (x) = 0, u (t) G dip (x (t)) x (0) = xo, v (0) = vo 
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where the first equality holds for almost all t G [0,oo[, and the inclusion holds for all t € 

Remark 4.5. The results obtained in this paper still hold if = Vif the gradient of a convex, continuously 
differentiable function fj : H —> K is replaced by a maximal monotone cocoercive operator B : H —> H 
(see [T] ). 


5 . Perspective 

Let us list some interesting questions to be examined in the future: 

(1) Study the asymptotic stability properties, corresponding to the case where T = +oo, in connection 
with the convergence results of [2] and [7]. 

(2) Extend the results to the case where the Levenberg-Marquart regularization term is given in a 
closed-loop form, X{t) = a(||i;(t)||) as in [6]. 

(3) Study the same questions for the corresponding forward-backward algorithms, see [T]. 
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